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Abstract 

We present a perturbative calculation of finite-size effects near T c of 
the ip lattice model in a d-dimensional cubic geometry of size L with 
periodic boundary conditions for d > 4. The structural differences 
between the ip 4 lattice theory and the p 4 field theory found previ- 
ously in the spherical limit are shown to exist also for a finite number 
of components of the order parameter. The two-variable finite-size 
scaling functions of the field theory are nonuniversal whereas those 
of the lattice theory are independent of the nonuniversal model pa- 
rameters. One-loop results for finite-size scaling functions are derived. 
Their structure disagrees with the single- variable scaling form of the 
lowest-mode approximation for any finite £/L where £ is the bulk cor- 
relation length. At T c , the large-L behavior becomes lowest-mode like 
for the lattice model but not for the field-theoretic model. Character- 
istic temperatures close to T c of the lattice model, such as T max (L) 
of the maximum of the susceptibility are found to scale asymp- 
totically as T c — T max (L) ~ L~ d / 2 , in agreement with previous Monte 
Carlo (MC) data for the five-dimensional Ising model. We also predict 
Xmax ~ L d l 2 asymptotically. On a quantitative level, the asymptotic 
amplitudes of this large -L behavior close to T c have not been observed 
in previous MC simulations at d = 5 because of nonnegligible finite- 
size terms ~ L^~ d ^ 2 caused by the inhomogeneous modes. These 
terms identify the possible origin of a significant discrepancy between 
the lowest-mode approximation and previous MC data. MC data of 
larger systems would be desirable for testing the magnitude of the 
£(4— d )/ 2 anc [ J_^~ d terms predicted by our theory. 



PACS numbers: 05.70.Jk, 64.60.-i, 75.40.Mg 
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1. Introduction 

A detailed understanding of the range of applicability of the ip 4 field theory in 
d dimensions is of fundamental interest to statistical and elementary particle 
physics jrj. The perturbative treatment of the critical behavior of the ip 4 field 
theory in d < 4 dimensions is known to be nontrivial because of the problem 
of infrared divergences. This problem is solved by the renormalization-group 
theory [l], §]. Above four dimensions where the critical behavior is mean- 
field like, no infrared problems of perturbation theory arise and no necessity 
exists for invoking the renormalization group. Thus the (p 4 theory above four 
dimensions appears to be free of essential problems. 

This is true, however, only for infinite systems. For the ip 4 field theory 
of confined systems above four dimensions [l], U there are several aspects of 
general interest — such as the nature of the fundamental reference lengths, 
the range of validity of universal finite-size scaling, the relevance of inho- 
mogeneous fluctuations, and the significance of lattice effects — that have 
remained unresolved until recently. These issues have turned out [4-6] to 
be closely related to the longstanding problem regarding the verification of 
earlier phenomenological |7|, |] and analytical predictions for d > 4 and 
regarding the various attempts to test these predictions by means of Monte 
Carlo (MC) simulations for the five-dimensional Ising model [7-12]. Clar- 
ifying these issues is of substantial interest for a better understanding of 
finite-size effects and of the concept of finite-size scaling [13-18], not only for 
d > 4 but also for the limit d — > 4. 

Recently [4-6] we have shown, on the basis of exact results in the limit 
n —>■ oo of the 0(n) symmetric if 4 theory, that finite-size effects for d > 4, 
for cubic geometry and periodic boundary conditions, are more complicated 
and less universal than predicted previously and that therefore the previous 
analyses of MC data were not conclusive. In particular we have found that 
lattice effects and inhomogeneous fluctuations of the order parameter play an 
unexpectedly import role and that two reference lengths must be employed 
in a finite-size scaling description. So far, however, no direct justification was 
given for our conclusions to be valid also for the more relevant case of lattice 
systems with a finite number n of components of the order parameter. 

It is the purpose of the present paper to provide this justification. We 
shall present a perturbative treatment of a ip 4 lattice model in one-loop or- 
der that leads to quantitative predictions of asymptotic finite-size effects for 
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d > 4 and n = 1. We shall show that the previous arguments M demonstrat- 
ing the necessity of using two scaling variables (rather than a single scaling 
variable) remain valid also for finite n for both the field-theoretic and the 
lattice model. We also confirm that the finite-size effects of the ip A lattice 
model differ fundamentally from those of the ip 4 field theory for general n. 
This implies that the Landau-Ginzburg- Wilson continuum Hamiltonian for 
an n component order parameter does not correctly describe the finite-size 
effects of spin models on lattices with periodic boundary conditions above 
the upper critical dimension. 

More specifically, we study the case of cubic geometry (volume L d ) with 
periodic boundary conditions and calculate the asymptotic finite-size scaling 
form of the order-parameter distribution function P($) where <£> is the spatial 
average of the fluctuating local order parameter ip. From P($) we derive 
the asymptotic finite-size scaling functions of the susceptibilty, of the order 
parameter and of the Binder cumulant ||19|| . Two scaling variables 

x = t(L/£ ) 2 , t = (T-T c )/T c (1) 

and 

y = {L/kf~ d (2) 

are needed where £o and Iq are reference lengths related to the bulk correlation 
length £, similar to the case n — > oo [4-6]. £o is the amplitude of £ for T > T c 
at vanishing external field h whereas 1$ is (proportional to) the amplitude 
of £ at T = T c for small h 0. The second length 1$ is connected with the 
four-point coupling u according to / ~ ul^ d ~ 4 \ As an alternative choice of 
scaling variables we also employ w and y where 

w = xy-W = t{L/i) d ' 2 (3) 

with the reference length 

I = Uio/loY /d ■ (4) 

In addition to the lengths £o, and L there is the microscopic length a or 
A , i.e., the lattice spacing of the lattice model or the inverse cutoff of the 
field-theoretic model. For short-range interactions, £o and Iq are expected to 
be of 0(a) or 0(A _1 ). Our scaling functions presented in Section 4 are valid 
in the asymptotic range L > a, ( > a or LA ^> 1, £A 3> 1. 
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The role of the length Iq is twofold. Since the dangerous irrelevant charac- 
ter of Mo [20, 21 1 exists already at the mean-field level, the length Iq appears 
via i in the variable w ~ UQ l ^ 2 tL d l 2 already at the level of the lowest-mode 
approximation which takes only homogeneous fluctuations into account JJ. 
The second important role of Iq originates from uq being the coupling of the 
inhomogeneous higher modes. Here uq does not have a dangerous irrelevant 
character. In fact, these higher modes are relevant for d > 4 as has been 
demonstrated for n — > oo || and will be shown in the present paper to be 
valid for general n, contrary to opposite statements in the literature [1,3]. 
The length Iq sets the length scale of the finite-size effects arising from these 
modes. 

Both scaling variables x and y or w and y must be employed in general, 
i.e., at any finite value of £/L < oo in the entire asymptotic L~ x — plane 
(Fig.l) to provide a complete description of asymptotic finite-size effects of 
the (p A theory. Our description is consistent with the general scaling struc- 
ture in terms of tL 2 and uoL A ~ d proposed by Privman and Fisher but 
inconsistent with the reduced structure proposed by Binder et al. [7] and 
with the lowest-mode approximation of Brezin and Zinn- Justin [3] in terms 
of a single scaling variable tL d / 2 equivalent to w. We find that it is only 
the region between the curved dotted lines in Fig. 1 where the single- variable 
scaling forms of Refs. [3] and [7] are justifiable for the lattice model, but 
not for the field-theoretic model. The region between the curved lines cor- 
responds to the special case £/L —>■ oo in the limit L — ► oo and \t\ — > at 
finite w where the large - L behavior becomes lowest-mode like for the lattice 
model. For t = (w = 0) this was found previously [4-6] for the case of the 
susceptibility and of the Binder cumulant, as conjectured in Ref. p2 . As a 



consequence we shall show that characteristic temperatures, in the sense of 
"pseudocritical" temperatures |]|] such as T max (L) of the maxima of the sus- 
ceptibility or the "effective critical temperature" || T C (L) where the magneti- 
zation has its maximum slope, scale asymptotically as T c — T max (L) ~ L~ d / 2 
or T c — T C (L) ~ L~ d l 2 for the y? 4 lattice model, in agreement with previ- 
ous Monte Carlo (MC) data |§. Correspondingly we predict Xmax ~ L d l 2 
asymptotically for the lattice model. 

On a quantitative level, our theory predicts that the large -L behavior 
close to T c is strongly affected by nonnegligible finite-size terms ~ 
caused by the higher (inhomogeneous) modes, as demonstrated recently for 
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the Binder cumulant at T c 0. In the analysis of Ref. [9] the observed 
"slow approach to the scaling limit" was considered to be the most signifi- 
cant discrepancy between the lowest-mode prediction || and the MC data. 
Our theory identifies the terms ~ jj^- d )l 2 as the possible origin of this dis- 
crepancy. We also show that, for the same reason, the successful method of 



determining bulk T c via the intersection point of the Binder cumulant |Tj| is 
not accurately applicable to finite spin models of small size in d — 5 dimen- 
sions, as demonstrated in Ref. [6]. New MC simulations over a larger range 
of L would be desirable for testing the predictions of our theory regarding 
the magnitude of the L^~ d ^ 2 and L A ~ d terms. 

In Section II we derive some of the bulk properties of the lattice model 
for n = 1 above four dimensions in one-loop order. In particular we identify 
the amplitudes of the correlation length at h = for T > T c as well as at 
T = T c for h ^ 0. In Section III we derive the effective Hamiltonian and 
the order-parameter distribution function in one-loop order. Applications 
to the asymptotic (large L 3> a, small \t\ <C 1) finite-size scaling functions 
and predictions of the lattice model for d — 5 are presented and discussed 
in Section IV. Results for the field-theoretic model are briefly presented in 
Section V. We summarize the results of our paper in Section VI. 

2. Lattice model: Bulk properties for d > 4 

We consider a (p A lattice Hamiltonian H for the one-component variables (fi 
with — oo < ifi < oo on the lattice points x« of a simple-cubic lattice in a cube 
with volume V = L d and with periodic boundary conditions. We assume 



H = a d \ ££V + u (tf) 2 - M + £ ^M<Pi ~ <Pjf (5) 

[ i i,j ) 

where a is the lattice spacing. The couplings .1^ are dimensionless quantities. 
The variables (fj have the Fourier representation 

w = j; d E e ik ' x ^ k ■ (6) 

In terms of the Fourier components 

^ = ~a d Y. e ^ X ^ ( 7 ) 
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the Hamiltonian H reads 

H = L~ d \ [ r o + M0k fi-k - hifo 
k 1 



+ U L 3d Y (0k 0k') (0k" 0-k-k'-k") (8) 
kk'k" 

where 

k = ~ 2 [J(0) - J(k)] (9) 
cr 

with 

J(k) = £ ^e- ik (Xl - x ^ . (10) 

The summation 5Z k runs over discrete k vectors with components 
kj = 2nmj/L, rrij = 0, ±1, ±2, ■ • ■ , j = 1,2, ■■■,(! with a finite cutoff 
—A = —7i /a < kj < it /a = A. We assume a finite-range pair interaction 
such that its Fourier transform has the small k behavior 

J k = J k 2 + 0{k]k 2 ^ (11) 

with 

Jo = \{a/L) d £(^/a 2 )(x, _ x .f . (12) 

The complete information on thermodynamic properties is contained in 
the Gibbs free energy per unit volume (in units of /c#T) 

/ = -^hxjvcpeM-H) (13) 

where the symbol / T>(p is the usual abbreviation for the multiple integral over 
the real and imaginary parts of (the finite number of) the Fourier components 

y?k- 

Recently we have found [[J in the large-n limit that the two reference 
lengths of the finite-size scaling functions for d > 4 are determined by the 
two bulk correlation-length amplitudes £o at h = for T > T c and Iq at 
T = T c for small h. We shall see that this property remains valid, as far 
as £o is concerned, also for n — 1 in one- loop order. With regard to Iq, 
an additional n-dependent prefactor arises that is 1 in the large-n limit and 
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g-i/22-i/3 £ or n = i [ n one- loop order (see Eq. (28) below). For this purpose 
we need to identify the amplitudes of the bulk correlation length for d > 4. 
We also calculate the bulk amplitudes of the order parameter M b and of the 
susceptibilities xt an d Xb above and below T c as reference quantities for the 
corresponding finite-size effects. 

All of our calculations are carried out at finite cutoff A (finite lattice 
spacing a). First we derive the asymptotic form of the bulk susceptibility xi 
at h = above and below T c as well as at T = T c for small h. The bulk 
Gibbs free energy density is denoted by f b = lim /. In terms of the bulk 

L^oo 

order parameter 

M b = - lim df/dh (14) 

L^oo 

the bulk Helmholtz free energy density T b = f b + M b h reads in one-loop 
order 

r b (r , M b ) = ir M 6 2 + « M fe 4 + ^ J ln(r + Ylu Q M 2 b + J k ) (15) 

k 

where J k stands for (2n)~ d J d d k with \ki\ < A. Above T c for h = 0, the 
inverse bulk susceptibility (xt)^ 1 is 

(xtr 1 = {d 2 T b /dM 2 b ) Mb= , = r + I2u J (r + J^ 1 + O(u 2 ) (16) 

k 

which determines the critical value r 0c of r as 



r 0c = - I2u J Jk 1 + O(u 2 ) . (17) 

k 

Thus we rewrite (x^) 1 above T c in terms of r — r 0c as 



(Xb) 1 = ( r o-r 0c 

where 



12u / J k 



/ 



+ O (u 2 ) (18) 



- r 0c = a t, t = (T - T c )/T c . (19) 

Note that the integral in Eq. (18) exists only for d > 4 and only for finite 
cutoff. The spontaneous value M s of the bulk order parameter for h — > 
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below T c is determined by dTb/dMb = 0. This yields for d > 4 



M 2 = (^"'(roc-ro) 



1 + 24m / J k 

k 



(20) 



The inverse susceptibility (Xb ) 1 for /i below T c is in one-loop order 
d 2 T\ 



(Xb)- 1 = 



9MiJ Mb=Ms 



2 (r 0c - r ) 



- 12u J 



T~ 2 



+ O(«o) • (21) 



From the equation of state at T ( 



/i = 4w M& 



— 36«o J 



T~ 2 



(22) 



we obtain the h dependence of the inverse bulk susceptibility x c 1 a ^ as 



y-l = ^± = 3h 2 ^U U 



1 — 36w y "Ac 

k 



1/3 



+ O(ug) . (23) 



For T > T c , the bulk susceptibility at finite wave vector q 

~2d 

X6(q) = lim — £ < ^ <Pj > e-^-^ (24) 

hJ 

has the one-loop form (for both h — and h ^ 0) 

x b (q) _1 = Xb(O)- 1 + J q [i + o(«g)] • (25) 

Thus the square of the bulk correlation length for T > T c is given by 

e = Xb (0) [<9x6(k)-7«9k 2 ] k=o = j Q X6 ( ) + («g) . (26) 

Substituting Eqs. (18) and (23) into Eqs. (25) and (26) yields the asymptotic 
form for d > 4 

£ = Sot- 1/2 , i > , h = , (27) 



and 



f = 3- 1 / 2 2- 1 / 3 / (/i 2 / d+ V - 1 )- 1 / 6 ,t = 0,h^O 



with the lengths 



„-l/2 rl/2 

a j 



1 + 12m / J k 2 



1/2 



+ 0(« 2 ) 



(28) 



(29) 



and 



u Jq 



1 + 36w / J k 

k 



_n V(«i-4) 



(30) 



These lengths will appear also in the finite-size scaling functions in Sect. 
IV. We see that for d > 4 the fluctuations (that enter via the one-loop 
integrals) only modify the amplitudes but do not change the mean-field t 
and h dependence. The "dangerous" uo dependence of £ at T c , Eqs. (28) 
and (30), is clearly seen in l ~ Uq 4 ^ We note that both £ and Iq are cutoff 
dependent via / k J k 2 . In rewriting [1 - 12m / k J k 2 ]~ 1/2 as [1 + 12u / k J k 2 ] 1/2 
in £ (and similarly in Z ) we have been guided by the resummed forms of £o 
and l in the limit n oo (Eqs. (141) and (142) in Ref. 

Corresponding results can be derived for the continuum f 4 Hamiltonian 
(see Eq. (62) below) with periodic boundary conditions, similar to the case 
n — > oo studied previously [Q. This amounts essentially to replacing J k by 
k 2 in the equations given above. As far as bulk properties of xtiXbi £ an d 
Mb are concerned, only the nonuniversal amplitudes are modified but the 
t and h dependence remains identical for the field-theoretic and the lattice 
ip A model. For the finite system, however, lattice effects become significant 
as we shall see in the subsequent Sections. For the specific heat even the 
(finite) bulk value at T c turns out to be different for the field-theoretic and 
the lattice model, similar to the case n — > oo [|J. 



3. Order-parameter distribution function for d > 4 
3.1. General form in one-loop order 

A perturbation approach to finite-size effects of the (p lattice model for d > 4 
can be set up in a way similar to the previous finite-size perturbation theory 



9 



for d < 4 above and below T c [ Z3"f . We decompose 

= $ + a i 



(31) 



and shall derive an effective Hamiltonian H e ^ [^j, for the lowest (homo- 
geneous) mode 



a x - 



(32) 



by a perturbative treatment of the fluctuations of the higher modes 

a 3 = L~ d ^ ^ . (33) 

Correspondingly we write the lattice Hamiltonian, Eq. (5), in the form 

H = H {$) + H 1 + H 2 , (34) 
with the lowest-mode Hamiltonian 



H ($) = L d (-r $ 2 + u $ 4 - h$) 



the Gaussian part 



with 



H x = L d ~( r 0L + Jk)°k 5"-k 



f 0L = r + 12w M 2 , 



and the perturbation part 

H 2 = ~a d ^[6n ($ 2 -M 2 )a| + 4 ? i a| 



(35) 



(36) 



(37) 



(38) 



Unlike the case d < 4, we must work, for d > 4, at finite cutoff. Thus we 
shall incorporate the finite shift r 0c ~ O(u ), Eq. (17), of the parameter 
r whereas in the previous [23] dimensional regularization at infinite cutoff 
there was no O(u ) contribution to r 0c . Thus we define 



M 2 



(39) 
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where now 



H, 



eff 



and 



= L a 



Z 



J d$exp(-# e// ) 



(40) 



(41) 



contain the shifted variable r — r 0c . 

The partition function is decomposed as 



(XI 

Z = JVipe~ H = J d$ exp[-(iJo + r)] 



(42) 



where 



r($) = -In | XVexp [-(//! + if 2 )] (43) 
is determined by the higher modes within perturbation theory. We rewrite 

h ($) + r($) = r(o) + # e// ($) (44) 

and define the order-parameter distribution function 

1 



P(<&) 



cxp 



Z*ff 

oo 

Z eff = J d$exp [-H eff ($) 



(45) 



(46) 



In a perturbation calculation with respect to H 2 we obtain the effective 
Hamiltonian in one-loop order 



H eff {$) 



r e ff $ 2 + u e ff $ 4 + 0($ 6 ) -m 



(47) 



where 



ro ff = r - r 0c + I2u 



-2 
k 



+lUulM*L- d J2(r 0L + J k )- 2 , 
k^o 



(48) 
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u e ff = u -36u 2 L- d J2(r 0L + -k)~ 2 ■ (49) 

In Eq. (48) we have added and subtracted r 0c and have replaced f 0L , in the 
0{uq) terms, by 

roL = r - r 0c + 12w M 2 . (50) 

This quantity is a positive function of r — r 0c for arbitrary — oo < r — r 0c < oo 
at any finite L. 

Moments of the distribution function can now be calculated as 

oo 

< $ m > = J rf$ $ m P($) (51) 



and 



oo 

< i$i m > = J rf$ i$i m p($) . (52) 



Because of the (one-loop) $ 4 structure of if 6 -^, these averages can be ex- 
pressed in terms of the well-known functions 

oo 

Jds s m exp(-iy S 2 - s 4 ) 

tmOO = ^ss (53) 

/ ds exp(-iFs 2 - s 4 ) 



that appear also in the finite-size theory below four dimensions p3| ^4 
The moments determine several thermodynamic quantities such as [23[ the 
susceptibilities 



X 



L d < $ 2 > , (54) 



X~ = L d {< $ 2 > - < |$| > 2 ) , (55) 
the "magnetization" M = < |$| > , and the Binder cumulant 

U = 1 - - < $ 4 > / < $ 2 > 2 . (56) 
3 
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In terms of the effective parameters rjj and u e J^ they can be expressed in 
one-loop order as 



X 



{L d /ul fS ) 1 ' 2 A (F e// ) 



(57) 



X 



(L d /u 



§ 2 (Y eff ) - $ t {Y efS f 



(58) 



M 



(L d ul Sf )- l ' A a (r e// ) 



(59) 



u 



with the dimensionless quantity 



yeff = L d/2 r e//(^// ) -l/2_ 



(60) 



(61) 



We note that at this stage of perturbation theory these expressions do not yet 
represent a systematic expansion with respect to the coupling u , compare 
Eqs. (5.19)-(5.22) of Ref. g|. 

Corresponding formulas are obtained for the case of the Landau-Ginzburg- 
Wilson continuum Hamiltonian 



H 



d d x 



V 



^r ^ 2 + -^(w) 2 + *W 4 - hip 



(62) 



with the field (p(x) by the replacement Jk — > k 2 in the sums of the one-loop 
terms of the effective parameters r$ and u e ^ . 

A justification of the above perturbation theory can be given in terms of 
the order-parameter distribution function where all higher modes are treated 



in a nonperturbative way [24] 



3.2 Asymptotic form of the effective parameters 

In order to study the asymptotic finite-size critical behavior we shall consider 
the limit of L > a, ^ > a or LA ^> 1, £A 3> 1. For this purpose we 



13 



decompose the perturbation part of r e J^ and u e J^ into bulk integrals and 
finite-size contributions in the following way, 



0"o - r 0c ) < 1 - I2u / J k (r 0L + J k 



+ U4u 2 M 2 



r L + Jk) 2 ~ I Jk(r L + Jk 



-i 



eff 



with 



12«oAi (rwO - 144«X A 2(r 0i ), 

m - 36mq / (r L + Jk)~ 2 + 36ulA 2 (r L) 



L- d £ (r 0L + J k )" 



A m (r 0L ) = y (r L + Jk)' 



(63) 
(64) 

(65) 



In the lowest-mode approximation we would have simply = ro , u'q 
Uq. Up to this point, the determination of the effective Hamiltonian for the 
field-theoretic model, Eq. (62), is still parallel to that of the lattice model. 
The corresponding formulas are simply obtained by the replacement — > k 2 . 
The crucial difference, however, comes from the large-L behavior of A m , as 
shown recently for the special case m = 1 and r L = 00. For general r L 
we find, for the lattice model, the cutoff-independent large-L behavior 

A m (r 0L ) = Jv m Im{ro L Jv l L 2 ) L 2m ~ d + 0{e- L l~ a ) 



with 



I m (x) = (2tt) 



-2m 



dy y 



m-l e -(xy/4TT 2 ) 



d/2 



K(y) d + 1 



(66) 



(67) 



where K(y) = J2 ex P(~yj 2 )- F° r the field-theoretic model, however, the 

j=-oo 

corresponding large-L behavior differs significantly according to the cutoff 
dependent result 



f(r 0L + k 2 )- m L- d J2(r OL + k 2 )- m = I m (roLL 2 )L 2m - d 

+ A d - 2m {a m (d,r 0L A- 2 )(AL)- 2 + 0[(AL)- 4 ]} (68) 
14 



where 
a m (d,r 0L A~ 2 ) 

This leads to 

„*ff 



d 



3(2n] 



d-2 



dye 



-y 2 x 



L-l 



d-l 



exp 



(l+r 0L A 2 )x 



(69) 



(r ~ r 0c ) U - 12txq / k 2 (r 0L + k : 



+ 144 M 2 M 2 ^ / (r 0L + k 2 



- 12w n 



k 2 (r 0i + k 2 



-i 



- 144w 2 M 2 



J 1 (r 0L L 2 )L 2 - d + A d - 2 a 1 Kr 0L A- 2 )(AL)- 2 ; 

I 2 (r 0L L 2 )L*- d + A d ^a 2 (d,r 0L A- 2 )(AL)- 2 } , (70) 



U JJ 



u - 36mq y (r 0L + k 2 ) 



2\-2 



+ 36m 2 



J 2 (r 0L L 2 )L 4 - d + A d - 4 a 2 (ci,r 0L A- 2 )(AL)- 2 (71) 



, d-4 



for the field-theoretic model. We note that for d < 4 the cutoff dependent 
terms in Eqs. (70) and (71) vanish in the limit A — > oo. For d > 4, however, 
part of these terms are divergent for A — >• oo and cannot be dropped. In par- 
ticular these terms carry the important size dependence ~ L~ 2 of r e J^ which 
is not present in Eq.(63) for the lattice model and which has been overlooked 
previously [0,0]. Employing the method of dimensional regularization in 
Eq. (68) would mean that the finite results for d < 4 at A = oo are continued 
analytically to d > 4. This would yield the same (cutoff-independent) result 
for r e J^ and u e J^ of the field-theoretic case as of the lattice model. Thus di- 
mensional regularization would omit the important analytic L~ 2 dependence 
in Eqs. (68) and (70). This omission, however, cannot be justified — unlike 
the omission of an analytic t dependence in bulk critical phenomena. Thus 
the method of dimensional regularization may yield misleading results in the 
finite-size field theory above the upper critical dimension. 

The remaining bulk integrals in r e J^ and u e Q have finite limits for tol — > 
(large L, small \t\) for both the lattice and field-theoretic model. Taking the 
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limit 7"ol — ► in these integrals is justified only if \t\ <C 1 and (L/a) _d//2 <C 1 
or (AL)~ d / 2 <C 1. This restriction should be kept in mind when applying our 
asymptotic scaling functions to MC data of spin models of small size. The 
asymptotic expressions for the lattice model for d > 4 read 



r eff 



1 - 12m 



-12u J^L 2 -%(r 0L J^L 2 ) 
-U4ulM*J^L A - d I 2 (r 0L J^L 2 ). 



eff 



U 



1 - 36u / J k 2 + 36u J^h{r QL J^L 2 )L A - d 



(72) 
(73) 



The corresponding results of the field-theoretic model for d > 4 are ob- 
tained from Eqs. (70) and (71) as 



eff 



{r ~ r 0c ) 



12u k" 



lUu 2 M 2 L 4 - d I 2 (r 0L L 2 ) 

,2 A rf-4 



-12«o L^hir^L') 



U4u 2 A d - 4 M 2 a 2 (d, r 0L A" 2 )(AL)- 



(74) 



WO 



1 - 36m / k" 4 + 3Gu L 4 ~ d I 2 (r 0L L 2 ) 
Jk 

+36M A d - 4 a 2 (d, r 0L A- 2 ) (AL) " 



(75) 



Substituting Eqs. (72) — (75) into Eqs. (47) - (49) completes our calcula- 
tion of the asymptotic form of H e ^ and of the order-parameter distribution 
function -P($), Eq. (45), in one-loop order for d > 4 and n — 1. The restric- 
tion "asymptotic" means that these results, Eqs. (72)-(75), are applicable 
to arbitrary r 0L L 2 only in the sense that L/d must be large and \t\ must be 
small. 

As found already in the large-n limit the leading "shift of T c " [3] 



2-rf 



in To , Eq. (72), is not just a temperature independent constant ~ L 
for the lattice model but a more complicated function of tqlL 2 . For the 
field-theoretic model the leading shift is proportional to L~ 2 according to 
Eq. (74) and is nonuniversal, i.e., explicitly cutoff-dependent. This result 
differs from the (temperature independent) shift ~ L 2 ~ d predicted for the 
field-theoretic model || and from the shift ~ L~ d l 2 considered in previous 
work [7—9]. Our shifts are caused by the (inhomogeneous) higher modes of 
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the order-parameter fluctuations. They cannot be neglected even for large 
L/a (except for the extreme case of the bulk limit) and cannot be regarded 
only as "corrections" to the lowest-mode approximation since they lead to 
a two-variable finite-size scaling structure for both the field-theoretic and 
the lattice model, in contrast to the single-variable scaling structure of the 
lowest-mode approximation, as will be further discussed in Section 4.2. 

Our results can be generalized to n > 1 by means of a nonperturbative 
treatment of the order-parameter distribution function of the 0(n) symmet- 
ric f 4 theory [24]. It is obvious that this does not change the conclusions 



regarding the structural differences between the finite-size effects of the field- 
theoretic and lattice versions of the (f A model. 

4. Finite-size scaling functions of the lattice model 

4.1 Analytic results 

In the following we consider only the case h = 0. Inspection of the asymptotic 
expressions of ro and Uq^ , Eqs. (72) and (73), shows that they depend on 
three different lengths £o, Eq.(29), Zq, Eq.(30), and L. Therefore there exist 
different ways of writing H e ^ in a finite-size scaling form. 

Considering the ratio £/L as a fundamental dimensionless variable [13-15] 
leads to the following scaling form 

# e // = F(z,x,y) = K eff {x,y)z 2 + u eff {x,y)z 4 , (76) 

where x and y are the scaling variables given in Eqs. (1) and (2) and z is 
the scaled order-parameter variable 

z = J " 1/2 L^ 2 ® . (77) 

The scaled effective parameters are in one-loop order 

r^f(x,y) = r e ff L 2 J^ = x - 12I x {f)y - UAd 2 {y Q )I 2 {f)y z l 2 (78) 

and 

u eff (x,y) = uf L'- d J^ 2 = y + 36I 2 (f)y 2 (79) 



with 



r 



z+120 2 (j/ o )2/ 1/2 (80) 
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and 

2/o = xy~ 

This leads to the finite-size scaling form 



1/2 



X 



L 2 PHx,y) 



and 



with the two-variable scaling functions 



U(x,y) 



Jo 1 



u 



effl 



-1/2 



x,y) ti 2 (Y(x,y)), 



(81) 

(82) 
(83) 

(84) 



J, 1 [u e "{x,y)Y 1/2 {MY(x,y)) - [#i{Y{x, v)))*} , (85) 

(86) 

(87) 
(88) 



J, 



-1/2 



eff (x,y) ~ 1/4 MY(x,y)), 



where 



l _l MY(x,y)) 



3[MY(x,y))} 



Y(x,y) 



~efft 



x,y) 



effl 



x,y) 



-1/2 



The traditional finite-size scaling theories for d < 4 [13-18] have no asymp- 
totic dependence on a second scaling variable y. We see that our scaling 
functions do not depend on the nonuniversal model parameters a, Jij, do, Uq, 
except via the length scales £o and Iq contained in x and y, and apart from 
the metric prefactors Jq 1 and Jq 1 ^ 2 in Eqs. (84)- (86). Thus we may consider 
these functions to be universal in a restricted sense, i.e., for a certain class 
of lattice models (rather than continuum models, see below). 

From the previous one-loop finite-size theory |23| and the successful com- 



parison with high-precision MC data in three dimensions |25j it has become 
clear that careful consideration must be devoted to the appropriate form of 
evaluating these one-loop results. The previous analysis indicated that the 



prefactor 



u 



eff 'W/2 



in Eqs. (84)- (88) should be further expanded with re- 



spect to the coupling uq, in the spirit of a systematic perturbation approach, 
see Eqs. (5.45), (5.46) and footnote 1 of Sec. 7 of Ref. p3| , as well as Eqs. 
(6.15), (6.16), (6.31) and (6.32) of Ref. |23|. Thus the expanded forms 



l u 



e//)-l/2 



y- 1/2 [l + lM 2 {f)y\- 1 or 



u 



eff\-l/2 



y- 1/2 [l - 18/ 2 (f)y] should 
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be substituted into Eqs. (84), (85) and (88), respectively [and similarly for 
(u e ff)~ l l A in Eq. (86)]. These expanded forms should be taken into account 
in a future quantitative comparison of Eqs. (84)-(88) with MC data. In the 
present paper we confine ourselves, for simplicity, to the unexpanded form of 
Eqs. (84)-(88), as has been done in the result for U(x,y) presented in Ref. 
[§]. The same comment applies to Eqs. (99)-(103) below. 

At T c (x = 0) we obtain from Eqs. (78) - (88) the large -L behavior in 
one-loop order for d > 4 

M c ~ L^/Zj-VZy-V^^^L-i/* , (90) 



lim U(Q,y) = 1 - -tf 4 (0)/tf 2 (0) 2 = 0.2705 . (91) 

L^oo 3 

The exponents in Eqs. (89) and (90) and the asymptotic value in Eq. (91) 
are identical with those obtained in the lowest-mode approximation at T c |J. 
The dangerous irrelevant character of uq ~ y is clearly exhibited in Eqs. (89) 
and (90) in the form of y~ x ^ 2 and y _1//4 . 

Alternatively we may employ, instead of x and y, the variables w and y 
where w is given by Eq. (3). This implies the following scaling form 

H eff = p(s,w,y) = ^r eff (w,y)s 2 + u eff (w,y)s 4 (92) 
with the scaled order-parameter variable 

a = Jo 1/2 lf~ A)lA $ . (93) 
The scaled effective parameters are 

r eff (w,y) = r eff y- 1/2 

= 12/i(f)y 1/2 - lU-8 2 (w)I 2 (f)y (94) 

and 

u e ff(w,y) = u eff y- 1 = l + 36I 2 (f)y (95) 
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with 



r = [w + \2$ 2 {w)]y 1/2 . 
This leads to the finite-size scaling forms 

x ± = L *l*p±fay) 

and 

M = L- d / A P M (w,y) 
with the two-variable scaling functions 

Pt(w,y) = A\u°"(w,y)]~ 1/2 2 (Y(u>,y)), 



(96) 

(97) 
(98) 

(99) 



p x (w, y) 
Puiw^y) 



A 



u e "(w, y)} ~ 1/2 {& 2 (Y(w, V)) - [MY(w, y))] 2 } , (100) 



u eff (w,y) ~ 1/A MY(w,y)) 



U(w,y) = l-±MY(w,y))/[MY(w,y)f 



(101) 
(102) 



where A = J l 



l,(4-d)/2 



and 



Y(w,y) 



r e ^(w,y) u eI:t (w,y) 



r, e fft 



-1/2 



(103) 



In the lowest-mode approximation the y dependence in Eqs. (92) - (103) is 
dropped and Eqs. (94) and (103) are replaced by Y = f e ^ = u 1 ^ 2 a tL d / 2 . 

These results will be discussed and applied to d — 5 in the following 
Subsections. 



4.2 Discussion 

Both sets of scaling variables (x,y) and (w,y) are useful in the analysis of 
the finite-size scaling structure. First we consider the (x, y) representation. 
In order to elucidate the effect of the fluctuations of the (inhomogeneous) 
higher modes above and below T c we assume y to be small (large L/l ) and 
expand x + an d X~ with respect to y at finite \x\ > 0, i.e., T ^ T c . This 
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yields 



X 



+ 



xt {l -^[x- 1 - h{x)] V - + 0{y 2 /x 2 )Y 



x>Q (104) 



X 



Xb {l + [iSx- 1 + 12/i(-2x) - 36x/ 2 (-2x) 



(105) 



where xt an d Xb are the bulk quantities given in Eqs. (18) and (21). Similar 
expressions can be derived for M and U. 

The terms ~ x~ l in the square brackets can be traced back to the lowest- 
mode contributions whereas the terms ~ and ~ I m (— 2x) arise from the 
higher (inhomogeneous) modes. If the latter terms were ignored one could 
rewrite m a lowest-mode form with the single variable x/y 1 ^ 2 , as noted 
previously in the case n — > oo [4—6]. For any finite \x\ = L 2 /!; 2 , however, i.e., 
along the straight dashed lines in the L^ 1 — plane (Fig. 1), there exists no 
argument that would allow one to ignore the I m terms arising from the higher 
modes. This proves the necessity of including two separate scaling variables 
in general. In particular it is misleading to consider the finite-size effects 
of the higher modes as a "correction" to the lowest-mode approximation 
— in the same sense as changes of mean-field exponents caused by critical 
fluctuations for d < 4 should not be considered as " corrections" . The crucial 
point is that the higher modes cause a new structure of the finite-size scaling 
functions that cannot be written in terms of a single variable x/y 1 ^ 2 (except 
for the special — > 0, y — > at finite x/y 1 ^ 2 , see below). This structural 
aspect is a matter of principle, regardless of how large or small the effect of 
the higher modes might be. 

In this context we take up a nontrivial aspect in the discussion in the 
previous literature about role played by the "shift of T". It was asserted 
that a term of the type ~ L 2 ~ d in the parameter r e J^ of H e ^ || represents a 
"correction to scaling" [9] or a "subdominant term" [12] that can be neglected 
in the large -L limit compared to the lowest mode part r — r 0c = ao^- This 
assertion is incorrect, however, for the reasons just given in the preceding 
paragraph. The term ~ L 2 ~ d in r e J^ , Eq. (72), is, in fact, the origin of 
the terms ~ I\{x) of x^ m Eqs. (104) and (105) which we have shown to 
represent nonnegligible contributions rather than "corrections". Similarly, 
the term ~ MqL 4- ^ in r$ , Eq. (72), is the origin of the nonnegligible 
higher-mode contribution 2x) to x~ m Eq. (105). 
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The expansion in Eqs. (104) and (105) breaks down in the limit of small 
\x\, i.e., large £/L. This includes the large -L limit at T = T c where the 
exponent of the susceptibility Xc ~ L d ^ 2 and the Binder cumulant U c have 
been found [4—6] to agree with the lowest-mode approximation for the lattice 
model (but not for the field-theoretic model). Here this result is seen from 
the representation in terms of w and y as given in Eqs. (92)- (103). In 
this representation the single-variable lowest-mode like structure appears as 
the leading w dependence whereas the higher-mode contributions ~ I\ and 
~ I2 are multiplied by y 1 ^ 2 and y. As noted previously @] these higher- 
mode contributions are not of a dangerous irrelevant character even though 
the dangerous irrelevant four-point coupling uq determines the length scale 



We shall see below that the y 1 ' 2 terms are quantitatively important. Nev- 
ertheless, at first sight it seems justified to consider the latter contributions 
as asymptotically negligible in the limit y — > corresponding to L — > 00. In 
the terminology of the renormalization group, this limit corresponds to ap- 
proaching the " Gaussian fixed point" of the dimesionless four-point coupling 
uoL 4 ~ d . Neglecting the ^/-dependence in this limit is justified, however, only 
if w is kept finite, i.e., if \t\ vanishes sufficiently strongly. Keeping w finite for 
L — > 00 is a special case corresponding to paths in the L~ x — plane where 
the ratio £/L diverges as L^ -4 ^ 4 . Such paths become asymptotically parallel 
to the vertical axis (Fig. 1). This includes the special case T = T c , L — > 00. 
The description of the entire L^ 1 — plane, on the other hand, requires 
both w and y as generic scaling variables in order to correctly include the 
w — > 00 limit which corresponds to a finite ratio £/L for L — ► 00. The latter 
limit w —>■ 00 cannot be taken correctly within the single-variable scaling 
structure such as P^(w,0) and within the reduced scaling form F(s,w,0). 
This structure does not capture the complete finite-size effects presented in 
Eq. (104) and (105) above. For this reason the inhomogeneous fluctuations 
must be considered as relevant, for finite L, in the sense of the renormal- 
ization group. For the same reason the reduced scaling form [equivalent to 
F(s,w,0)} that was proposed by Binder et al. [0 for general £/L (not only 
for £/L = 00), is not valid. 

Our order-parameter distribution function P($) ~ exp(—H e ^) can be 
compared with the zero-field probability distribution function of Binder || 
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T9| below T c 



Pl{s) 



const \ exp 



- [8 



m b YL d /2 Xb 



+ exp 



-(s + m b yL d /2 Xb (106) 



where m b ~ l^l 1 ^ 2 and Xb ~ are L independent bulk quantities. In Pl(s) 
the temperature dependence enters in the form [L/l(t)] with the "thermody- 
namic length" [8,9,19] l(t) ~ Itl -2 ^. Our theory identifies the relevant length 
scale i of the corresponding variable w, Eq. (3), in terms of a combination 
of l and £o according to Eq. (4). The distribution function Pl{s) has been 
invoked as an argument in support of the single-variable scaling structure of 
the free energy (at h — 0) proposed by Binder et al. J7|. For finite £/L, 
however, our results do not agree with the structure of Pl(s) which does 
not contain the important y dependence reflected in the shift ~ I\l?~ d in 
Eqs. (72), (78), (94) that is caused by the inhomogeneous modes. Thus 
the double Gaussian form of -Pl(s), Eq. (106), as well as the underlying 
theory of Gaussian thermodynamic fluctuations [p6[| , are not applicable to 
finite systems with periodic boundary conditions in the critical region above 
the upper critical dimension. This is remarkable in view of the fact that 
mean-field theory becomes exact in the bulk limit for d > 4. 



4.3 Predictions for d = 5 

We illustrate and further discuss our results for the lattice model for the case 
d = 5 which can be compared with MC data of the five-dimensional Ising 
model. In Fig. 2 we plot the order-parameter distribution function in terms 
of F, Eqs. (76)-(81), 

P($, t, L, u ) d$ = ^M-F^x.y)] ^ ^ (iQ7) 

/ dzexp [-F(z,x,y)] 

— oo 

for typical values of x and y above, at and below T c . The shape of these 
functions resembles that of the MC data shown in Fig. 1 of Ref. ||. In 
order to demonstrate the effect of the higher modes on these functions we 
show the order parameter distribution function in Fig. 3 in terms of F, Eqs. 



23 



(92)-(96), 



P($, t, L, u ) rf$ = ^eM-Hs^y)} ^ m 

J dsexp —F(s,w,y) 



with given w but for several values of L/l Q including the limiting function 
for L/l — > oo at fixed w. The latter function has the structure Po(&) of 
the lowest-mode approximation (with one-loop expressions for the reference 
lengths £ and l ). This does not mean, however, that Po{&) is the exact 
representation of P in the large -L limit in general. The constraint w = 
const < oo restricts the validity of Po($) only to the special region \t\L d l 2 < 
oo corresponding to £/L — > oo (region between the curved dotted lines in Fig. 
1). The width of this region in the L~ l — £ _1 plane vanishes asymptotically for 
L — > oo. This special region is of interest because it contains characteristic 
("pseudocritical" [13]) temperatures close to T c such as T max (L) and T C (L) 
to be defined below. 

The scaling functions i~, Eqs. (99) and (100), of \ + an d X~ correspond- 
ing to the order-parameter distribution function of Fig. 3 are shown in Fig. 
4. Similar plots can be made for M and U. For comparison with MC data we 
refer to Figs. 11-14 of Ref. [9]. At first sight, the changes due to the variation 
of L/Iq appear to be small. At the level of accuracy of previous MC data [|], 
however, these changes and their disagreements with the lowest-mode pre- 
dictions || have been clearly detected and have been considered as a major 
discrepancy, regardless of their smallness, because of their unexplained weak 
L dependence. 

Here we further elucidate the deviations from the lowest-mode predictions 
by plotting in Fig. 5 the scaling functions of U, x ± and M at and below T c as 
functions of the reduced length L/Iq. As originally found in Ref. for the 
example of the Binder cumulant at T c , the slow approach to the asymptotic 
(L — > oo) values arises from the y 1 ^ 2 ~ L( 4_c 0/ 2 terms. This slow approach 
was observed in the previous MC data and, at that time, gave sufficient 
reason to doubt the correctness of the lowest-mode predictions || which do 
not have a weak subleading L-dependence at T c . Our theory now shows 
that subsequent attempts [10-12] to explain the discrepancies did not resolve 
the problems. In particular, the bulk form of the renormalization-group flow 
equations employed by Blote and Luijten Jl0| led to an apparent confirmation 



of the (incorrect) shift ~ L 2 d predicted in Ref. j3| for the field-theoretic 
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model. These bulk flow equations |I0| do not correctly describe the finite-size 
effects of the ip A lattice model either. We believe that our theory identifies the 
origin of the previous discrepancy, apart from possible quantitative aspects 
which we shall address elsewhere, after a quantitative identification of the 
lengths £ and l . 

An interesting consequence of the existence of the limiting function Po(<$>) 
mentioned above is the existence of limiting scaling functions such as P~ (w, 0) 
of x~ for L/Iq — > oo at fixed w [Fig. 4 (b)]. For finite L, x~ exhibits a maxi- 
mum Xmax below T c at a temperature T max (L). The asymptotic L dependence 
of T c —T max (L) can be inferred from the fact that P^(w, 0) has a temperature 
dependence only of the form of w ~ tL d l 2 . This implies the large -L behavior 
T c — T max (L) ~ L~ d l 2 and correspondingly Xmax ~ L d ^ 2 . Similar arguments 
lead to our prediction T C — T C (L) ~ L~ d l 2 where T C (L) is the "effective critical 
temperature" || at which the magnetization has its maximum slope. The 
same power law is valid for the temperature at which the specific heat has its 
maximum. These power laws ~ L~ d l 2 agree with the MC data ||. The true 
asymptotic amplitudes of these power laws, however, have not been observed 
in previous MC simulations in d = 5 dimensions because of the slow approach 
of the subleading terms ~ 2 7 ( 4 ~ a! )/ 2 towards L — ► oo mentioned above. MC 
simulations of larger systems would be desirable for testing the magnitude of 
such subleading terms predicted by our theory. 

As indicated already in Figs. 3 and 4 of Ref. [6], we also point to 
an additional interesting effect of practical importance. In Fig. 6 we have 
plotted the Binder cumulant as a function of x for several values of L/Iq. 
Without the effect of the slowly decaying contribution ~ y 1 ^ 2 ~ £( 4 ~ d )/ 2 one 
would have expected || a well identifiable intersection point of these curves 
if L is, say, larger than 10 a. For d < 4, this features has been a standard 
and successful empirical method of determining the value of bulk T c from MC 
data of finite systems. Our previous |J and present figures demonstrate that 
this method is not accurately applicable, without additional information, to 
systems with d — 5. 



5. Finite-size scaling functions of the ip field theory 

For the field-theoretic model the effective parameters rp and u e ^ , even 
in their asymptotic form in Eqs. (74) and (75), depend explicitly on the 
length A -1 , in addition to the lengths £o,^o and L, as found already in the 
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limit n — > oo M. This means that none of the original nonuniversal model 
parameters a , uq and A becomes unimportant even close to T c and for large 
L. In a scaled form the effective parameters read 

r eff = r eff L 2 = x _ u _ 144 ^ 2 (y )/ 2 (f)y 3 / 2 

- 12u A d - 4 a^d^A- 2 ) 

144( Mo A d - 4 ) 3/ ¥ 2 (yo)a 2 (rf,rA- 2 )(AL)- d / 2 (109) 

and 

+ 36( Mo A d - 4 ) 2 a 2 (d,fA- 2 )(AL) 2 - d (110) 

where f and ?/o are given in Eqs. (80) and (81). The last term ~ L 2 ~ d in Eq. 
(110) can be neglected asymptotically. 

Substituting these expressions into Eqs. (76), (84) — (88) and (107) 
(with J = 1) yields the finite-size scaling functions of the order-parameter 
distribution function and of the quantities x^^M and U. The two- variable 
finite-size scaling functions depend on x and y and, in addition, explicitly on 
the nonuniversal parameter uoA d_4 . Thus the scaling functions are nonuni- 
versal for n=l, and obviously also for general n. At T c , the asymptotic 
power laws of % an d M are found to be 

xt(L) = L 2 P+(0,y)~L d - 2 , (111) 

M C (L) = LM/ 2 P M (0,y) ~ L' 1 (112) 

for the field-theoretic model which differ from those of the lattice model in 
Eqs. (89) and (90). The asymptotic value of U at T c for the field-theoretic 
model is in one-loop order 

lim U(0,y) = 2/3 (113) 

which is far from that of the lattice model in Eq. (91). The significant 
differences between Eqs. (89)-(91) and Eqs. (111)-(113) are due to the L- 
independent but cutoff-dependent term ~ u A d ~ 4 in Eq. (109), similar to 
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the constant additive term in Eq. (122) of Ref. M. Existing MC data for 
Ising models on d = 5 lattices (such as the MC result Umc = 0.319 ± 0.017 
in Ref. |5|) clearly disagree with these field-theoretic results and rule out 
the possibility that the ip 4 field theory provides a correct description of finite 
lattice systems above the upper critical dimension. In particular, the predic- 
tion of a breakdown of universality for finite systems above the upper critical 
dimension constitutes a serious failure of the continuum approximation for 
lattice systems. The present results for finite n confirm our earlier [4-6,22] 
assertion regarding the applicability of the ip 4 field theory for d > 4. 



6. Summary and conclusions 

We summarize and further comment on the results of this paper as follows. 

On the basis of a one-loop calculation for the ip 4 model on a lattice and 
for the ip 4 continuum model in a cubic geometry with periodic boundary 
conditions above four dimensions we have shown that our general conclusions 
regarding universality and finite-size scaling inferred from the large-n limit [4- 
6] remain valid for finite n. In particular, ip 4 field theory based on the Landau- 
Ginzburg- Wilson continuum Hamiltonian |1[ does not correctly describe the 
leading finite-size effects of spin systems on a lattice with d > 4. 

Although the critical exponents of mean-field theory are exact for bulk 
systems above four dimensions, the thermodynamic theory of Gaussian fluc- 
tuations PBj is not applicable to finite systems with periodic boundary con- 



ditions in the critical region for d > 4. 

Finite-size scaling in terms of a single scaling variable, as predicted by 
the phenomenological theory of Binder et al. and by the lowest-mode 
approximation of Brezin and Zinn- Justin ||, is not valid for the ip 4 field 
theory. For the p> 4 lattice model it is not valid for any finite £/L where £ is the 
bulk correlation length. As originally conjectured in Ref. [p2]| , lowest-mode 
like large -L behavior is asymptotically correct for the lattice model at T c as 
shown previously for the susceptibility [4—6] and the Binder cumulant |J; 
furthermore, it is valid in the small region of finite \w\ ~ \t\L d l 2 in the large 
-L and small \t\ limit (Fig. 1). This region, corresponding to a divergent ratio 
£/L ~ £( d ~ 4 )/ 4 — » oo, represents only a small part (between the two curved 
dotted lines of Fig.l) of the general finite-size scaling regime (of arbitrary 
finite £/L) for which earlier theories |3], [7| were originally thought to be 
valid. Our two-variable finite-size scaling structure is consistent with that 
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proposed by Privman and Fisher [^1|] but is significantly less universal than 
anticipated previously [15]. 

The inhomogeneous higher modes have been shown to be relevant above 
the upper critical dimension, contrary to different statements in the previous 
literature [1,3,10,12,15,27-43]. We have identified the characteristic length 
scale Iq, Eq. (30), of the finite-size effects of the higher modes in terms of the 
amplitude of the bulk correlation length £ at T = T c for small external field 
h. The one-loop finite-size effects arising from the relevant higher modes do 
not represent " corrections" to the lowest-mode approximation but constitute 
a generic part of the correct finite-size scaling structure. By contrast, two- 
loop contributions are expected to represent only quantitative corrections 
that will not change the scaling structure. 

The " shift of T c " || in the temperature variable rp of the exponential 
order-parameter distribution function is proportional to L~ 2 for the field- 
theoretical model and proportional to I(t, L~ 1 )L 2 ~ d for the lattice model 
where the function I(t,L~ l ) has a finite limit 7(0,0). The effects caused by 
these shifts remain nonnegligible at any finite ratio £/L even in the large 
-L limit as demonstrated in Eqs. (104) and (105) for the example of the 
susceptibility above and below T c . 

The " shift of T c " ~ L 2 ~ d mentioned in the preceding paragraph must be 
distinguished from shifts of characteristic temperatures, in the sense of pseud- 



ocritical temperatures [n|, such as the temperature T max (L) at which the sus- 
ceptibility x~(t,L) has its maximum, or the "effective critical temperature" 
T C (L) where the order parameter has its maximum slope. We find that 
these "shifts" have the asymptotic (large L) behavior T c — T max (L) ~ L~ d l 2 
and T C — T C (L) ~ L~ d l 2 for the (p 4 lattice model. Similarly our theory implies 
Xmax ~ L d / 2 asymptotically. This is a simple consequence of the fact that 
the order-parameter distribution function shown in Fig. 3 has a finite limit 
for L —>■ oo at finite w and that the position of T max (L) and T C (L) remain 
located in the temperature region of finite w in the limit L — > oo. 

Our theory identifies the possible origin of a significant discrepancy be- 
tween MC data at d = 5 [9] and the lowest-mode prediction for the Binder 
cumulant at T c |3j in terms of slowly decaying finite-size terms ~ 
(Fig. 5). These terms also mask the true asymptotic amplitudes of the 
power laws Xc ~ L d l 2 and M c ~ L~ d l A at T c . For the same reason the 
method of determining bulk T c (from MC data via the intersection point of 
the Binder cumulant) is demonstrated in Fig. 6 to become quantitatively 
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inaccurate at d = 5, as found originally in Ref. ||. 

Quantitative predictions for various asymptotic finite-size scaling func- 
tions have been made for d = 5 and n = 1 (Figs. 2—6). These predictions 
are expected to be valid for sufficiently large L/a and small \ t\. The true range 
of applicability remains to be explored by quantitative comparisons with MC 
data, after an appropriate identification of the nonuniversal lengths £o an d 
Iq. As noted previously H, it is not yet established whether the ip 4 model on 
a finite lattice is fully equivalent to finite spin models regarding the leading 
and subleading finite-size effects. 

Our results indicate that in the limit d — 4 — > + , for systems with 
periodic boundary conditions, different amplitudes of finite-size effects at 
d = 4 are obtained depending on whether a lattice model or a continuum 
model is considered. In view of this possible ambiguity at d = 4 the limiting 
behaviour for 4 — d — > + (i.e., e — > + in the standard e = 4 — d expansion ) 
should also be reexamined for lattice models at finite lattice spacing and for 
continuum models at finite cutoff. 

Note added 

After completion of the present work we received a preprint " Finite- 
size scaling above the upper critical dimension revisited: The case of the 
five-dimensional Ising model" by E. Luijten, K. Binder, and H.W.J. Blote 
where the authors compare the asymptotic result for U(x,y) in the (unex- 
panded) form of Eqs. (87) and (88) with their Monte Carlo data of the five- 
dimensional Ising model. The authors confirm the " occurrence of spurious 
cumulant intersections" predicted in Figs. 3 and 4 of Ref. and agree with 
the slow convergence of finite-size effects for L —>■ oo found in Ref. |J which 
essentially resolves the longstanding discrepancies noted in the MC studies 
in Refs. [7-9] regarding the Binder cumulant. On a much more quantitative 
level than considered previously [3-12], however, the authors estimate the 
length l as Iq = 0.603 (13) and claim to find new "significant discrepancies" 
between their MC data for small L and our asymptotic (large L) one-loop 
result for \ + m the (unexpanded) form of Eq. (84). 

We doubt the significance of the quantitative deviations for small system 
sizes L = 4 and 8 shown in their Figs. 7(a) and (b), except for the sign and 
curvature of the deviations from the large-L behavior of the susceptibility 
shown in their Figs. 8 and 9. We propose that the latter issue can be 
resolved essentially on the basis of our complete non- asymptotic one-loop 



29 



expression for if e " presented in Eqs. (47)- (49) of this paper or on the basis 



of the underlying order-parameter distribution function [24|, rather than by 



an asymptotic two-loop calculation suggested by Luijten et al. We note that 
similar non-asymptotic effects are well known for small spin systems at T c in 
three dimensions as discussed in the context of Fig. 14 of Ref. [18]. 

We doubt the reliability of the estimate of l = 0.603 (13) by Luijten et 
al. since it was found by applying the asymptotic (L — ► oo) expression for 
xt in their Eq. (31) to non- asymptotic MC data. Our Fig. 5b indicates that 
the apparent (6 percent) mismatch between theory and MC data for L < 22 
in Fig. 9 of Luijten et al. is due to this inadequate estimate of 1$. 

Part of the remarks by Luijten et al. regarding the limiting case t — > 0, 
L — > oo at fixed tL d l 2 agree with our earlier and present independent findings. 
We disagree, however, with their claim that "there is no contradiction at all" 
between our finite-size theory and the ideas of Ref. J7J. First, we note that 
the ideas of Ref. [7] fail for the continuum ip A model. Second, we maintain 
that the single- variable scaling structure (for h = 0) proposed in Ref. does 
not capture the correct structure of finite-size effects of the lattice model at 
any finite value of £/L (see Fig.l). In particular the ideas of Ref. [7] do not 
lead to a correct description of the finite-size departures from bulk critical 
behavior at small but fixed \t\ in the asymptotic range < \t\ < 1 [compare 
Eqs. (104) and (105)] to which an acceptable finite-size scaling structure 
should be applicable (such as the general structure proposed in Ref. pl|). 
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Figure Captions 

Fig.l. Asymptotic L^ 1 — plane (schematic plot) above (t > 0) and 
below (t < 0) T c for the lattice model above four dimensions where L is 
the system size and £ is the bulk correlation length (in units of the lattice 
spacing). The straight dashed lines correspond to paths at constant finite 
ratio £/L. The curved dotted lines represent paths at constant finite tL d ^ 2 for 
L — ► oo and t — > corresponding to a divergent ratio £/L ~ £( d_4 )/ 4 — ► oo. 
The single-variable scaling form in terms of w = t(L/l) d l 2 , Eqs. (3) and (4), 
is valid only in the region between the curved dotted lines. The two-variable 
scaling form in terms of x = tL 2 /£q and y = {L/Iq) or w and y is necessary 
in the entire L" 1 — plane where £/L is finite. 

Fig. 2. Theoretical prediction of the scaled order-parameter distribution 
function P(3>, t, L, Uq)L'^ 2 ^ 2 J^ 2 of the v? 4 lattice model in ci = 5 dimen- 
sions vs z, Eq. (77), in the form of Eq. (107) where F(z,x,y) is given by 
Eqs. (76)— (81). The parameter values are x = 0.25 above T c (a), x = at 
T c (b), and x = —0.75 below T c (c), where \x\ = L 2 /£ 2 . The reduced length 
is y^ 1 = L/Iq— 16. Compare Fig. 1 of Ref. 0. 

Fig. 3. Theoretical prediction of the scaled order-parameter distribution 
function P($, t, L, uo)L~ d ^ 4 J^ 2 Zq 4 of the y? 4 lattice model in d = 5 
dimensions vs s, Eq. (93), in the form of Eq. (108) where F(s,w,y) is given 
by Eqs. (92) — (96). The parameter values are w = 1 above T c (a), w = at 
T c (b), and w = -3 below T c (c) where w = t{L/l) d ' 2 . The curves are shown 
for two reduced lengths L/l = 8 (dotted lines) and L/l = 32 (dashed lines) 
as well as for the limiting case L/l = oo (solid lines) at fixed w. 

Fig. 4. Theoretical prediction of the finite-size scaling functions 
P+(w,y)Jol ( Q 4 ^ 2 and P~(w,y)JolQ d 4 ^ 2 of the susceptibilies x + ; Eq. (99), 

(a) and x~i Eq. (100), (b) of the lattice model in d = 5 dimensions vs w = 
t(L/l) d l 2 for several values of the reduced length L/l , corresponding to Fig. 
3, including the limiting case L/l = oo at fixed w. In this representation the 
position of the maximum of P~ attains a finite value for L/Iq — > oo which 
determines T c — T max ~ L~ d l 2 . Compare Fig. 13 of Ref. [9]. 

Fig. 5. Theoretical prediction of various finite-size scaling functions 
of the lattice model in d = 5 dimensions, U(0,y) (a) , P£(0,y)J o lQ 

(b) , P M {Q,y)jl /2 lt m (c), P-(-3,y)J lt 4)/2 (d) , P"(-6, y) Jo/S d " 4)/2 (e) 
according to Eqs. (99)-(103), as a function of the scaled length y~ x = L/l . 
The arrows indicate the asymptotic one-loop values for L —>■ oo. Compare 
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Fig.2 of Ref.[6]. 

Fig. 6. Theoretical prediction of the Binder cumulant U(x,y), Eq. (87), 
as a function of the scaled reduced temperature x = t(L/£ ) 2 , Eq. (1), in 
the range —4 < x < 4 (a) and in the range —0.4 < x < 0.4 (b), for several 
values of the scaled length y~ x = (L// ) d ~ 4 , at d = 5 : L/l = 8 (dotted 
line), L/l Q = 16 (dot-dashed line), L/l = 32 (dashed line). The bulk value 
[thin solid line in (a)] is 2/3 below T c and above T c . The cross indicates 
the asymptotic one-loop value U(0, 0) = 0.2705 at T c , L — > oo, Eq. (91). 
Compare Figs. 3 and 4 of Ref. |§. 
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